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We consider plane waves propagating in quadratic nonlinear slab waveguides with nonlinear quasi- 
phase- matching gratings. We predict analytically and verify numerically the complete gain spectrum 
for transverse modulational instability, including hitherto undescribed higher order gain bands. 
OCIS codes: 190.4410, 190.4420 



I. INTRODUCTION 

With the maturing of the quasi-phase-matching 
(QPM) technique, in particular by clcctric-ficld poling 
of ferro-electric materials, such as LiNb03 , by poling 
of polymers 13] and quantum-well disordering in semi- 
conductors [3j , the number of applications of quadratic 
nonlinear (or x^) materials has increased significantly. 
Even complicated QPM grating structures are now com- 
mercially available in periodically poled LiNbC>3 (PPLN). 
It is therefore more important than ever to have com- 
plete knowledge of the effects a QPM grating has on the 
properties of x' 2 ' materials. The most fundamental ef- 
fect of a QPM grating, with a certain spectrum of spa- 
tial wave vectors, is to allow noncritical phase matching 
at all wavelengths for which the wave-vector mismatch 
Afc matches a component of the grating spectrum. Thus 
QPM gratings allow for efficient multiple wavelength sec- 
ond harmonic generation (SHG) 0], which may, for ex- 
ample, be used for multiple-channel wavelength conver- 
sion p. 

In addition to providing phase matching, QPM grat- 
ings have the basic effect that they induce an asymmetric 
cubic nonlincarity (ACN) in the equations for the aver- 
age field in the form of self- and cross-phase modulation 
terms §. This ACN appears in linear and/or nonlinear 
periodic QPM gratings of arbitrary shape pj , it can be fo- 
cusing or defocusing, depending on the sign of the phase 
mismatch 0, and its strength can be increased (e.g., 
dominating the Kerr nonlincarity) by modulation of the 
grating 8]. In continuous- wave operation the ACN in- 
duces an intensity-dependent phase mismatch, just as the 
inherent Kerr nonlinearity, with potential use in switch- 
ing applications 0, ^j. The ACN further explains a 
wide spectrum of novel fundamental properties of soli- 
tons [1 0E1 and modulational instability (MI) [rl[T^|. 
Importantly the ACN is a general effect of non-phase- 
matched wave interaction and as such appear also in 
homogeneous x^ materials in the cascading limit. In 
fact, in this case the asymmetric signature of the ACN 
may be measured as the difference between the proper- 
ties in upconversion and downconversion, since there is 
no effective x' 2 ' nonlinearity competing with the ACN 
as in QPM gratings. Such an experiment was recently 



reported 0] and thus the ACN has now been confirmed 
both numerically and experimentally. 

In this paper, we present the results of a complete 
study into the modulational instability of beams in 1+1D 
QPM waveguides. We find that in general, the MI gain 
spectrum has a multipeaked structure of up to three fun- 
damental bands with accompanying overtones. Our pre- 
vious work[H G3 concentrated on the structure of the 
fundamental (long-wave instability) bands at low trans- 
verse wavenumbers. These bands are due to MI of the 
averaged field, and are predicted by averaged equations, 
provided the effective ACNs are taken into account for 
the most accurate description. For example, for QPM 
with a simultaneous linear and nonlinear grating and/or 
with a nonlinear grating with a dc- value (as with QPM in 
polymers and semiconductors 0), the ACNs can sup- 
press these fundamental bands, making plane waves mod- 
ulationally stable over hundreds of diffraction lengths 
[l2L IT^ | . Exact Floquet calculations and numerical sim- 
ulations confirmed the predictions for the low-frequency 
bands, but also revealed high-frequency bands not pre- 
dicted by the averaged theory. These gain bands were 
surmised to be related to the inherent instability in homo- 
geneous (non-phase-matched) x' 2 ' media [l2L fl3|. How- 
ever, an accurate analytic description of the first overtone 
bands was not derived and the higher overtone bands 
were not discussed at all. 

Here we concentrate on PPLN slab waveguides, which 
are x^ media with a purely nonlinear QPM grating with 
no dc-component. Although more-general kinds of grat- 
ing can be analysed with our method, we choose to fo- 
cus on the simpler case for clarity. For example, with 
this simple QPM grating we do not need to take into ac- 
count the ACN for an accurate description of MI [l2[ Il3| , 
however our analytic method does allow for such cubic 
nonlinearities to be included. We present the first an- 
alytical and numerical description of the complete MI 
spectrum with all overtones. We show explicitly that the 
overtone series are caused by MI in the rapidly varying 
components of the propagating fields, which are in turn 
induced by the grating. We derive approximate analytic 
expressions for the positions of the gain bands and com- 
pare them with exact Floquet results and direct numeri- 
cal simulation, to find good agreement. 
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II. METHOD 

We consider a linearly polarized electric field propagat- 
ing in a lossless x^ slab waveguide with a QPM grating 
under conditions for (the most efficient) first-order QPM 
and type I SHG. The normalized dynamical equations for 
the slowly varying envelopes of the fundamental E±(x,z) 
and second harmonic (SH) E2(x,z) take the form 



dE 1 ld 2 E 1 , s , , N 

i^- + -^r^-+xi.z)ElE 2 e^(i(3z) = 0, (1) 



dz 2 dx 2 

.dE 2 ld 2 E 2 
j 1 

dz 4 dx 2 



+ X (z)EfeM-if3z) = 0, (2) 



where the asterisk means complex conjugate. The x and 
z scales have been normalised by Xq and zo=kiX 2 l , respec- 
tively, and f3—(k 2 — 2k±)zo is the normalized wave- vector 
mismatch. The nonlinearity is periodic in z with expan- 



d n exp(innz), 



(3) 



where d n —d*_ n (x is real) and the grating wave number 
k is defined to be positive. The grating will force the 
same periodicity in the propagating fields. We therefore 
expand these also: 



E\(x, z) — w n (x, z) exp (inKz), (4) 

n 

E 2 (x,z) — ^""^ v n (x, z) exp (inKz — i(3z), (5) 
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FIG. 1: Left: Maximum gain versus o n and e n for |6 n | 2 =l 
and |c| 2 =2. Right: a n (y n ) for v n real and found from Eq. 1181 
(solid line), from Eq. 1191 with plus sign (dotted line), and 
Eq. 1191 with minus sign (dashed line), for /3=0, A=l, and 
k=670. 



w n (x,z) = [w n + e n {x, z)] exp (iAz), (8) 



v n (x,z) = [v n + fJLn(x, z)] exp («2Az). (9) 



Substitution into Eqs. I|6I7|) gives the linearized equations 



(L[ - nK)e n + dn+m,-l-»{w%Jii + vi£* m ) = 0,(10) 

l,m 

(L' 2 - TlKjfln + 2 X! d n-m-l+sW m €i = 0,(11) 



where Lj = Lj — jA. Writing the perturbations in the 
form 



where [3=[3 — sk is the residual mismatch, with s=sign(/3) 
for first-order QPM. Substituting all expansions into 
Eqs. (11I2H gives 



e n (x,z) — S^(z) exp (ivx) + 5^*(z) exp (— ivx%12) 
fjt„(x, z) = 6^'(z) exp (ivx) + 5^*(z) exp (— ivx)(13) 



(Li - nn)w n + ^2 d n+m -i- s w* m vi = 0, (6) 
(L 2 - nK)v n + ^ d n - m -i +s w m wi = 0, (7) 

where Lj — id z + d 2 /(2j) + (j — This set of equa- 
tions has plane- wave solutions of the form w n (x,z) = 
w n exp (iAz), v n (x, z) = v n exp (2iAz) dill. 

To study MI we consider small perturbations of the 
plane- wave solutions: 



one obtains a linear matrix equation for the perturba- 
tion 5 n =(Sn\ Sn\ 6n\ 5^) T , which couples all Fourier 
components. 

To derive a simple result we consider nearly phase- 
matched interaction with |/c|~|/3|3>l in the typical (for 
PPLN) square QPM grating, for which d n =0 for n even 
and d n =2/(iirn) for n odd. The requirement |k| ~ 
|/3 3> 1 allows us to assume that the amplitudes w n 
and v n vary slowly compared to exp (ikz) and that the 
higher harmonics are much smaller than the dc-fields, 
\u>o\, l^o I ^ |w„^o|: l^n^ol- 111 this case the structure of 
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and e n reveals three gain bands, with extrema remaining 
close to a n = — e„ (diagonal branch DB) and a n e n = 2|c| 2 
(hyperbolic branches HB + , HB_), as illustrated in Fig.Q] 
for \b n \ 2 = 1 and |c| 2 = 2. The extrema a n — —e n of DB 
correspond to 



v 2 = -8n K /3 - 4A + 4/3/3, 



(18) 



while the extrema a n e n = 2|c| 2 of HB + and HB_ corre- 
spond to 



v z n = 2/3 - 3hk - 5A ± J (2/9 - jik - 3A) 2 + 16|c| 2 . (19) 



FIG. 2: MI gain for A=l (dashed) and A=-l (solid) cal- 
culated by Floquet theory. Crosses and diamonds show the 
theoretical prediction 1151 . n—p=670. 

the matrix equation means that the evolution of the per- 
turbation in each component is to a good approximation 
decoupled from the other components, 
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6 n = M n 5 n , (14) 



where a n =—v 2 /2 — A — nn, b n =d,2 n -sVo, c—d- s WQ, and 
e n =— v 2 /A + f3 — 2 A — tik. The eigenvalues of M n are 



A 2 = A n ± v^ 2 - B n , (15) 

where 

^ - ^(l&n| 2 -a'-4|c| 2 - e 2 ), (16) 

Bn = (a„e„-2|c| 2 ) 2 -e 2 |&„| 2 . (17) 

Any positive real part of an eigenvalue corresponds to MI 
with the gain g n (y) = 5R(A„). 

III. RESULTS 

For the gratings that we consider here, |c| 2 ~ A(2A— (3) 
and \b n \ 2 ~ A 2 /(2n — s) 2 . Analysing the gain versus a n 



For n^O the HB + bands appear at v 2 ~ — Ann — 8A + 
4/3 and the HB_ bands at v 2 ~ — 2nn — 2A. Thus we have 
the structure of up to 3 gain bands in the average field 
(n=0), each with a set of equally spaced (in v 2 ) weak 
overtone gain bands (n ^ 0). However, the gain is large 
only near a„=0 (see Fig. 0. For n ^ 0, a n ~ can 
be satisfied without violating the assumptions only on 
HB_, and thus the overtones of the other branches will 
be small. 

Consider a system at exact phase-matching (3 = 0, with 
dimensionless parameters |A| = 1 and n = 670. For a 
LiNbOs QPM grating with period A g = 17.6/im designed 
to phase-match at wavelengths of around A = 1550nm, 
these dimensionless parameters correspond to having an 
input intensity of J = 6.7GW/m 2 , if one assumes a 
waveguide depth of yo = 3.1^m. The distance scaling pa- 
rameters for this case are xq = 21.6/mi and zq = 1.89mm. 
For fixed distance scaling, different values of the magni- 
tude of A correspond to different input intensities, and 
different signs of A correspond to different ratios of the 
fundamental to second harmonic [13J. On the other hand, 
keeping |A| fixed while varying the input intensity corre- 
sponds to changing the distance scaling: to a good ap- 
proximation, zq oc Z _ 2 , xq oc and hence rea/J. 

For A=l our analytic results predict only one funda- 
mental n=0 gain band (HB + ) and three gain bands in 
all overtones. However the overtones of DB and HB + are 
too small to be seen. Thus the visible gain bands have 
maxima at ^o=1.18 (HB + ), i/_i=36.64, ^_2=51.84, etc. 
(all HB_). In physical units these are v$ = 0.055//m -1 , 
V-i = 1.697/im~ 1 , and v-2 = 2.401/J,m _1 , respectively. 
The anal ytic al results agree with exact Floquet calcula- 
tions 0, , as seen in Fig. [5] 

For A=— 1 two fundamental n=0 bands are predicted 
(DB and HB + ), with one set of visible overtones from 
HB_. The maxima are at u =2 (DB), t/ =3.38 (HB+), 
^_i=36.70, ^_2=51.88, etc. (all HB_). In physical units 
these are i/ = 0.093/im _1 , v Q = 0.0157/im -1 , = 
1.700/im -1 , and z/_2 = 2.403/im~ 1 , respectively. This is 
again confirmed by the exact Floquet calculations seen in 
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seeded by noise also confirmed our results (see Fig. 01. 

We see from Eqs. (118119(1 that the dimensionless spec- 
tral positions of the higher-order bands depends on k, 
whereas the fundamental bands do not. The different de- 
pendence implies a different scaling with input intensity, 
since for fixed |A|, both the distance scalings x , zq and 
the dimensionless wave vector k depend on the intensity. 
The result is that the spatial frequencies of the higher- 
harmonic instabilities remain the same, in physical units, 
whereas the frequencies of the fundamental-band insta- 
bilities generally increase with the quartic root of the 
intensity. 



IV. SUMMARY 



FIG. 3: Gain calculated by numerical simulation [(a) and 
dotted line in (b-d)] and by Floquet theory [solid line in (b- 
d)]. k = /3 = 670 and A = -1. 



Fig-El Figure|21also shows the gain profiles predicted by 
Eq. (|15fl . which agree well with the exact Floquet results. 
Direct numerical simulation El of Eqs. (1112ft with MI 



In conclusion, we have presented a simple theory able 
to accurately predict the complete MI spectrum in gen- 
eral QPM gratings in x materials. In particular, we 
have predicted and verified that overtone gain bands orig- 
inate from MI in the higher-order Fourier components of 
the field. This research is supported by the Danish Tech- 
nical Research Council (Grant No. 26-00-0355) and the 
Australian Research Council. 
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